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Abstract 

The hadronic production of a Higgs boson (H) in association with b jets will play an important 
role in investigating the Higgs-boson couplings to Standard Model particles during Run II of the 
CERN Large Hadron Collider, and could in particular reveal the presence of anomalies in the 
assumed hierarchy of Yukawa couplings to the third-generation quarks. A very high degree of 
accuracy in the theoretical description of this process is crucial to implement the rich physics 
program that could lead to either direct or indirect evidence of new physics from Higgs-boson 
measurements. Aiming for accuracy in the theoretical modeling of R -|- 6-jet production, we have 
interfaced the analytic Next-to-Leading-Order QCD calculation of Hbb production with parton- 
shower Monte Carlo event generators in the POWHEG BOX framework. In this paper we describe the 
most relevant aspects of the implementation and present results for the production of R -|- 16 jet, 
H -\-2b jets, and R with no tagged 6 jets, in the form of kinematic distributions of the Higgs boson, 
of the 6 jets, and of the non-6 jets, at the 13 TeV Large Hadron Collider. The corresponding code 
is part of the public release of the POWHEG BOX. 
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I. INTRODUCTION 


The production of a Higgs boson with b jets at the CERN Large Hadron Collider (LHC) 
can provide essential information on the Higgs-boson couplings to third generation quarks, 
in particular to the bottom quark. Indeed, all the leading parton-level production processes 
{gg —)■ Hbb, qq —?■ Hbb, and bg —?■ Hb) involve a Higgs boson radiated from an external 
bottom quark, while at the loop-level the Higgs boson can also originate from internal loops 
of both bottom (leading) and top quarks (see, e.g., the discussion in [HIS]). In the Standard 
Model, the Higgs-boson couplings to both fermions and gauge bosons are just proportional 
to the particles’ masses, causing the Higgs-boson associated production with bottom quarks 
to be largely suppressed with respect to the major production mechanisms, like gluon-gluon 
fusion (mediated by a loop of top quarks) or vector-boson fusion and associated production 
with vector bosons (where the Higgs boson couples to W or Z bosons). The production with 
b jets is then further suppressed by the identihcation cuts usually placed on b jets for tagging 
purposes. This scenario can however be drastically different if the hierarchy of Higgs-boson 
Yukawa couplings is modihed by factors that typically enter in models with extended Higgs 
sectors, like Two Higgs Doublet Models. In the quest for unveiling the origin of the breaking 
of the electroweak gauge symmetry, the evidence for (or absence of) H + &-jet production 
at Run H of the LHC can therefore provide an essential piece of the puzzle. 

In view of its crucial role for the Higgs-boson physics program of Run H of the LHC, H+b- 
jets production has received quite some attention in the context of the LHC Higgs Cross 
Section Working Group M, and both ATLAS and CMS have used the Hbb production 
channel in all major studies to constrain supersymmetric models and other extensions of the 
Standard Model [HHIO]. On the theoretical side, it is essential to control and improve the 
accuracy with which we can estimate rates for H production with one or two b jets, and, 
with this regard, quite some progress has been made in the last few years. In the following 
we will briefly summarize the main results obtained in this context, while we refer to the 
existing literature for more exhaustive explanations and details. 

As for all processes involving multiple scales (such as the masses of the bottom quark or 
the Higgs boson, and Mh, as well as scales determined by the kinematics of collisions at 
TeV center-of-mass energies), the QCD perturbative prediction of iL -|- 6-quark production 
can be affected by the presence, at all orders, of large corrections proportional to logarithms 
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of ratios of these mass scales (e. g., log(m;,/Q) where indicatively one can assume Q ~ Mh or 
larger). The occurrence of these enhanced logarithmic corrections depends on the signature 
studied as well as on the kinematic regime considered. It can be shown that these corrections 
can be reabsorbed in the perturbative dehnition of a bottom-quark parton distribution 
function, and from this observation originates the prescription to calculate processes like 
H + b quarks in a 5 Flavor Scheme (5FS) where the b quark is treated as a light parton 
and can appear in the initial state. At hxed perturbative order this is an alternative to the 
usual 4 Flavor Scheme (4FS) calculation where only four light-flavor parton densities are 
assumed while the bottom quark is treated as massive and can only appear in the hnal state. 
The production of H + lb jet can be induced at tree level by bg —)■ bH in the 5FS, and by 
QQi 99 Hbb in the 4FS; while the production of H + 2b jets can only arise at lowest order 
via qq, gg —>■ Hbb and is therefore an unambiguous 4FS prediction. The two approaches 
only correspond to a different reordering of the (same) perturbative expansion, and their 
predictions tend to agree better the higher the perturbative order, showing the expected 
well behaved convergence of QCD predictions. Several discussions of this issue can be found 
in the literature, where the general 4FS/5FS formalism is thoroughly analyzed [HI [I2] as 
well as specialized to the particular case of H + 6-quark production m El El [13]. 

Next-to-Leading Order (NLO) QCD corrections to H + b jet production have been cal¬ 
culated both in the 5FS [Tl| and in the 4FS [T51 - IT7] . The hrst set of corrections is nowadays 
part of the NLO QCD 5FS prediction of H + lbiet [Sl iTTlITS] . while the second set of correc¬ 
tions are included in the NLO QCD 4FS prediction of both H + lb jet im and H + 2b jets 
(smsiiiE]. NLO QCD hxed-order results for both Hb and Hbb hadronic production can also 
be obtained via any of the public NLO automated tools such as MadGraph5_aMC@NL0 [T9] . 
GoSam [20[, or OpenLoops [21] ■ Of course, the totally inclusive cross section (with no b 
tagging) can be calculated in either scheme, and dedicated studies which includes up to 
Next-to-Next-to-Leading-Order (NNLO) QCD corrections have been presented in the liter¬ 
ature [H El El ESI E2H21]. 

In order to improve the accuracy of theoretical predictions for total cross sections 
and distributions, these NLO hxed-order results need to be consistently interfaced with 
parton-shower Monte Carlo generators like PYTHIA [25l ES] and HERWIG [271 ES] using 
one of the methods proposed in the literature, namely MCSNLO [221 ED] and POWHEG [STL 
ISS] . as implemented in specihc frameworks like, e. g., MadGraph5_aMC@NL0 [TDl 1ST] , the 
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POWHEG BOX [35], and SHERPA [36|- The implementation of Higgs-boson production with b 
quarks in MadGraph5_aMC®NL0 has been discussed in Ref. [2], where total cross sections have 
been given for both inclusive and exclusive production and distributions have been shown in 
particular for the inclusive case (no fe-j^t tagging). In this paper we present the implemen¬ 
tation of H + 5-jet production in the POWHEG BOX, based on the 4FS NLO QCD calculation 
of Hbb hadronic production of Ref. [16]. While Ref. [2] considers both the 5FS and a 4FS 
cases, we will only consider the 4FS case since we aim at presenting in particular results for 
both H + lb jet and H + 2b jets in the same framework. We note that the implementation of 
5-initiated processes in a NLO QCD parton-shower Monte Carlo is still being studied and, to 
our knowledge, it is not routinely available in any of the aforementioned frameworks. Hence 
our decision to only implement the 4FS case. The details of the implementation will be 
presented in Section [TT| Results for H + 15-jet, iL-P25-jet, and H with no tagged 5 jet will be 
given in Section using a specihc setup, for the purpose of illustrating the kind of studies 
that are now possible within the POWHEG BOX framework. Our conclusion are presented in 
Section [IVl 


II. IMPLEMENTATION 

The implementation of the Hbb process in the framework of the POWHEG BOX can be per¬ 
formed along the same lines as the related Hti process that has been considered in Ref. 071. 
While the POWHEG BOX package provides all process-independent building blocks, it requires 
a list of all independent flavor structures for the tree-level contributions at Leading Or¬ 
der (LO) and NLO, the Born and real-emission amplitudes squared, the hnite parts of the 
virtual contributions, the color- and spin-correlated amplitudes squared, and a parametriza- 
tion of the phase space for the Born process. The flavor structures and tree-level amplitudes 
can most conveniently be generated with the help of a tool based on MadGraph 4 [38113^ 
that is provided in the POWHEG BOX. The virtual contributions for the pp —)■ Hbb process 
are extracted from the NLO-QCD calculation of [16] and adapted to the format required 
by the POWHEG BOX. All building blocks are implemented in the 4FS, i.e. no contributions 
from incoming bottom quarks are taken into account and the bottom-quark mass is always 
considered to be non-zero. 

While at LO and in the real-emission contributions only diagrams including a Hbb cou- 
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pling emerge, in the virtual corrections also loop diagrams with a Hti coupling contribute. 
These are fully taken into account in the representative results discussed in this work. The 
user of the POWHEG BOX implementation can choose to de-activate the contributions including 
a top-quark Yukawa coupling via a switch in the input hie. This allows to rescale separately 
the two contributions as necessary to calculate Hbb production in, for instance, supersym¬ 
metric extensions of the Standard Model, as discussed in detail in [1], where also a rescaling 
prescription is provided. 

We note that, contrary to the case of Hti production, where the heavy-quark mass is typ¬ 
ically renormalized in the on-shell renormalization scheme, in the case of Hbb production the 
renormalized bottom-quark mass is often dehned in the MS renormalization scheme mm- 
While both renormalization schemes are perturbatively equivalent at NLO with differences 
only due to higher-order contributions, Hbb production processes have been found to be 
quite sensitive to the renormalization scheme via the bottom-mass dependence of the over¬ 
all bottom-quark Yukawa coupling. Indeed, as higher-order corrections beyond the one-loop 
level are partly taken care of, physical observables are often found to exhibit a better per¬ 
turbative behavior when the MS scheme is used for the the bottom-quark Yukawa coupling. 
Taking this into account, the results presented below have been obtained in the MS scheme. 
When using the POWHEG BOX implementation of the Hbb process, however, the user is free 
to choose either the on-shell or the MS renormalization scheme for the bottom-quark mass 
that enters the Yukawa coupling by setting the respective parameter in the input hie. 

Although in principle not necessary for obtaining hnite results, technical cuts at the 
generation level can help to improve the performance of the Monte-Carlo integration. For 
the computation of observables with identihed b jets we therefore recommend the use of a 
small cut on the transverse-momentum of the bottom quarks (e.g. p™* = 0.1 GeV) when the 
phase-space integration is performed. We have checked that hnal results for the respective 


scenarios in Section III do not change when generation cuts of p™* = 0.1 GeV or p™* = 1 GeV 
are imposed compared to the case where no generation cuts are applied. 

In order to verify the POWHEG BOX implementation of the Hbb process, we have performed 
a detailed comparison of cross sections and distributions at LO and NLO as obtained in the 
POWHEG BOX with the hxed-order code of [16], and found full agreement for all considered 
observables. In addition we have successfully compared our results to those of Ref. |2|. 
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III. RESULTS 


The code we developed is available from the webpage of the POWHEG BOX project, 
http://powhegbox.inib.infn.it/, With this version of the code the user is free to study 
Hbb production at a hadron collider in a customized setup. Here, we wish to discuss represen¬ 
tative results for Hbb production at the LHC with a center-of-mass energy of y/s = 13 TeV. 
We use the four-flavor MSTW2008 set of parton distribution functions iQl SI] as imple¬ 
mented in the LHAPDF library ^2], with the associated value of a*, a Higgs-boson mass 
oi rriH = 125 GeV, and a top-quark mass of irtt = 173 GeV. The on-shell mass of the bot¬ 
tom quark is set to = 4.75 GeV, resulting in an MS mass of mhimf,) = 4.34 GeV at 
NLO QGD. 

For the renormalization (/ir) and factorization (pp) scales we consider two options: hrst 
we use a fixed scale. 


/^o 


mjj -|- 2mb 

4 


( 1 ) 


and second a dynamical scale. 


Fo = ^ -, (2) 

i 

where the summation runs over the masses and transverse momenta of the Higgs boson and 
the partons in the hnal state of the hxed-order calculation. To assess the scale dependence 
of our predictions, we vary the renormalization and factorization scales, pr = and 
/xp = ipo simultaneously in the range ^ = 1/2 to .^ = 2. 

We have matched the fixed-order NLO calculation with PYTHIA-6.4.25. In order to be 
able to focus our discussion of the NLO-PPYTHIA results on genuine parton-shower effects, 
we did not activate multi-parton interactions, underlying event effects, or decays of the 
Higgs boson in the Monte-Garlo program, although each of these effects could in principle 
be accounted for by setting the respective parameters in PYTHIA. 

In our numerical analysis we consider the two scenarios with a Higgs boson produced in 
association with one or two identihed b jets, as well as the case with no tagged b jets. Jets 
of any type are reconstructed with the anti-fcp algorithm as implemented in the FAST JET 
package |13], with R = 0.5. A jet that contains either a bottom quark or antiquark, or a H 
meson, is considered a b jet. For our H -|- 26-jet analysis, we require at least two b jets with 
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a minimum transverse momentum in the central region of pseudorapidity, 

jet) > 25 GeV , | p(6—jet)| < 2.5 , (3) 

while for the H + 16-jet analysis only one identihed b jet fulhlling the above criteria is 
required. We do not impose any cuts on extra jets, unless stated otherwise. 

In Figs. [Tjj^we illustrate the impact of the parton shower on the hxed-order NLO QCD 
results for the transverse-momentum (p^) and pseudorapidity (p) distributions of the Higgs 
boson and the hardest of the identihed h jets in H + 26-jet production, respectively, for both 
the hxed- and the dynamical-scale choice /ip = Pr = po of Eqs. 0 and ([^, respectively. In 
Figs. |3]j^we show the corresponding distributions for if-|- 16-jet production, and in Fig. [^we 
show the Pt{H) and p(if) distributions for the inclusive case, i.e. when no 6 jets are tagged. 
We hnd that parton-shower effects do not signihcantly change the hxed-order NLO results of 
these distributions within the given statistical uncertainty in most of the kinematic regimes 
shown. Only in the Higgs pt distributions for all signatures, ff -|- 0,1, 26 jets, we hnd that 
parton-shower ehects decrease (enhance) the NLO results for small (large) values of pt{H). 
These distributions also show that the ehects of the parton shower do not signihcantly diher 
for a hxed- and dynamical-scale choice within their respective statistical uncertainties. 

In Fig.(left), we show the ehect of the parton shower on correlations of the two identihed 
6 jets in the H -|-26-jet case, in particular their invariant mass distribution {M{bb)) and their 
separation in the azimuthal-angle-pseudorapidity plane (f2(6, 6)). The impact of the parton 
shower is again not signihcant within the statistical errors (apart for very small values of 
Mibb)). 

To illustrate the behavior of distributions including a non-6 jet we show in Fig. the 
i?(6, j) distribution in the H -|- 26-jet case, where f?(6, j) is the separation of the hardest 
6 jet and the hardest non-6 jet in the azimuthal-angle-pseudorapidity plane. Here, we only 
consider non-6 jets with a transverse momentum larger than 25 GeV in the rapidity region of 
the detector \yj\ < 4.5. In the l.h.s. plot we show the comparison between hxed-order NLO 
results and results obtained after the same calculation is interfaced with PYTHIA6 in the 
POWHEG BOX framework, for the central value of both the hxed and dynamical scale. Since 
the ehect of parton shower and in particular of scale dependence (hxed vs dynamical) seems 
much bigger than for other distributions, we further investigate the scale dependence of the 
distribution, which is shown in the r.h.s plot of Fig. Glearly, the /2(6,j) distribution 
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FIG. 1: The pT (left) and 7] (right) distributions of the Higgs boson in ff + 26-jet production 
at NLO-QCD with no parton shower (blue short-dashed line for fixed scale, black long-dash- 
dotted line for dynamical scale), and with parton shower as obtained through PDWHEG-I-PYTHIA6 
(red medium-dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The 
lower panels show the ratios: R = ficT(NLO)/(i(T(P0WHEG-|-PYTHIA) for a fixed (blue points) and a 
dynamical (black points) scale, respectively. The error bars indicate the statistical uncertainties of 
the Monte-Carlo integration. 
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is affected by a large scale uncertainty, both for a fixed- and a dynamical-scale choice. 
This is typical of observables that are indeed effectively described only at LO by a given 
NLO calculation. In the case of Hbb production, the hardest non-6 jet can only stem 
from the real-emission contributions of the hard matrix element (namely qq, gg —)■ Hbb -|- g 
or qg,qg —)■ Hbb + q/q), or from extra radiation due to the parton shower. Therefore, 
distributions of the hardest non-6 jet are effectively described only with LO accuracy and 
are affected by a typical LO scale uncertainty. Gaining full NLO control on jet distributions 
with smaller scale uncertainties would require an NLO calculation for pp —)■ Hbb + jet. The 
different behavior for hxed- and dynamical-scale choices which appears in the l.h.s. plot 
of Fig. (where only the central value of each scale is used) is then due to the large scale 
uncertainty encountered in distributions effectively only described to LO accuracy in the 
hxed-order calculation. We hnd large differences between the hxed-order predictions with 
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FIG. 2: The pT (left) and r] (right) distributions of the hardest identified b jet in + 26-jet produc¬ 
tion at NLO-QCD with no parton shower (blue short-dashed line for fixed scale, black long-dash- 
dotted line for dynamical scale), and with parton shower as obtained through PDWHEG-I-PYTHIA6 
(red medium-dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The 
lower panels show the ratios: R = ficT(NLO)/(i(T(P0WHEG-|-PYTHIA) for a fixed (blue points) and a 
dynamical (black points) scale, respectively. The error bars indicate the statistical uncertainties of 
the Monte-Carlo integration. 


^ = 2 and ^ = 0.5 when choosing Rr = Rf = with po defined in Eqs. 0 and ([^. These 
differences are reduced once the hxed-order calculation is combined with the parton shower 
in the P0WHEG-PPYTHIA6 predictions. 

Indeed, to illustrate the effect of the parton shower on non-6-jet observables we show in 
Fig. I^the transverse momentum of the hardest non-6 jet, but this time with no cut applied 
on the non-6 jet. If no pt cuts are imposed on the non-6 jets, in a hxed-order calculation 
the transverse-momentum distribution of a non-6 jet becomes very large at small values of 
Pt, which is due to large contributions from the emission of partons of very soft or collinear 
type. This behavior is tamed once the POWHEG-Sudakov factor is applied, as it is the case in 
the POWHEG+PYTHIA result shown in Fig.j^for the case of inclusive and H + 2h-]ei production. 
Since this effect is not sensitive to the tagging of b jets, similar shapes are encountered in 
both analysis scenarios (and also in the H + 16-jet scenario which is not explicitly shown 
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FIG. 3: The pT (left) and rj (right) distributions of the Higgs boson in H + 16-jet production 
at NLO-QCD with no parton shower (blue short-dashed line for fixed scale, black long-dash- 
dotted line for dynamical scale), and with parton shower as obtained through PDWHEG-I-PYTHIA6 
(red medium-dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The 
lower panels show the ratios: R = ficT(NLO)/(i(T(P0WHEG-|-PYTHIA) for a fixed (blue points) and a 
dynamical (black points) scale, respectively. The error bars indicate the statistical uncertainties of 
the Monte-Carlo integration. 


here). 

In order to assess the theoretical uncertainties associated with the choice of renormaliza¬ 
tion and factorization scale for NLO distributions, we have computed the pt and rj distri¬ 
butions of the Higgs boson for all three signatures, i.e. ih -|- 0,1, 2b jets, and the pt and rj 
distributions of the hardest identihed b jet for both H + 16-jet and H + 26-jet production, 
for different choices of scales as discussed earlier. The corresponding results, obtained us¬ 
ing our P0WHEG-I-PYTHIA6 implementation, are shown in Fig. (no 6 tagging). Figs. 10 11 


{H + 16 jet), and Figs. 12 13 {H -|- 26 jets). The scale dependence of the results is consider¬ 
able, amounting to about ±25% in most regions of phase space. Using a hxed scale rather 
than a dynamical scale helps in slightly reducing the scale uncertainty of the NLO±PYTHIA 
results in all observables, especially at larger values of Prib), Pt{H) and in the central pseu¬ 
dorapidity region, although the effect is moderate. We note that also the total cross sections 
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FIG. 4; The pT (left) and r] (right) distributions of the hardest identified b jet in H + 16-jet produc¬ 
tion at NLO-QCD with no parton shower (blue short-dashed line for fixed scale, black long-dash- 
dotted line for dynamical scale), and with parton shower as obtained through PDWHEG-I-PYTHIA6 
(red medium-dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The 
lower panels show the ratios: R = ficT(NLO)/(i(T(P0WHEG-|-PYTHIA) for a fixed (blue points) and a 
dynamical (black points) scale, respectively. The error bars indicate the statistical uncertainties of 
the Monte-Carlo integration. 


in the three scenarios considered here exhibit a large scale uncertainty, for instance we find 
cr = 0.477 pb ± 18% for the total inclusive cross section obtained with our NLO-I-PYTHIA6 
implementation in the setup of Fig. [^for a hxed-scale choice. 


IV. CONCLUSIONS 

In this article we have presented the implementation of the NLO QCD calculation for 
Hbb production at a hadron collider (from Ref. [16]) in the POWHEG BOX package. We em¬ 
phasize how having H -|- 6-jet production available in the POWHEG BOX provides a crucial 
element of consistency for experimental studies that rely on the same framework for a broad 
variety of signal and background processes. The code is made publicly available so that it 
can be used for further studies of Hbb production at the LHC in the SM and in extensions 
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FIG. 5: The pT (left) and rj (right) distributions of the Higgs boson in the inclusive case at NLO- 
QCD with no parton shower (blue short-dashed line for fixed scale, black long-dash-dotted line 
for dynamical scale), and with parton shower as obtained through P0WHEG-I-PYTHIA6 (red medium- 
dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The lower panels show 
the ratios: R = d(T(NLO)/ficT(POWHEG-|-PYTHIA) for a fixed (blue points) and a dynamical (black 
points) scale, respectively. The error bars indicate the statistical uncertainties of the Monte-Carlo 
integration. 


of the SM with modified Yukawa couplings of the third generation quarks. Here, we consid¬ 
ered the SM and presented numerical results at hxed perturbative order and at NLO QCD 
matched with PYTHIA for selected representative setups. In particular, we discussed theo¬ 
retical uncertainties due to factorization/renormalization scale choices and variations and 
parton-shower effects for two analysis scenarios with a Higgs boson in association with one 
or two b jets, respectively. A sample of results for the case of inclusive production have also 
been presented. We found that parton-shower effects do not give rise to large distortions of 
observables related to the Higgs boson or identified b jets in iL -|- 16-jet and H + 26-jet pro¬ 
duction processes at the LHC. As expected, more pronounced effects occur in distributions 
involving non-6 jets, as shown, for example, by the transverse-momentum distribution of 
the hardest non-6 jet or by the distribution of the separation (i?(6, j)) between the hardest 
6 jet and the hardest non-6 jet. We studied the impact of different scale choices on various 
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FIG. 6; The M{bb) (left) and R{b,b) (right) distributions in H + 26-jet production at NLO- 
QCD Tvith no parton sho'wer (blue short-dashed line for fixed scale, black long-dash-dotted line 
for dynamical scale), and with parton shower as obtained through P0WHEG-I-PYTHIA6 (red medium- 
dashed line for fixed scale, orange long-dash-dotted line for dynamical scale). The lower panels show 
the ratios: R = d(T(NLO)/ficT(POWHEG-|-PYTHIA) for a fixed (blue points) and a dynamical (black 
points) scale, respectively. The error bars indicate the statistical uncertainties of the Monte-Carlo 
integration. 


distributions and found that the associated theoretical uncertainties can be considerable. 
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